home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 005 / epistat.arc / MCNEMAR.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1983-08-18  |  1.4 KB  |  34 lines
  1. 1  REM                      MCNEMAR'S TEST
  2. 2  REM             Written by Tracy L. Gustafson, M.D.
  3. 3  REM            Round Rock, Texas. Version 2.0, 1983
  4. 5  DEF SEG=&H40
  5. 6  A=PEEK(&H17): IF NOT(A AND &H20) THEN POKE &H17,(A AND (NOT &H20)) OR &H20
  6. 7  DEF SEG: KEY OFF
  7. 10  CLEAR: DEFINT A-C,N,T,Z: SCREEN 0,0: WIDTH 80: COLOR 7,0,4
  8. 20  CLS: PRINT TAB(20);"KEY";STRING$(36,205);"CLOSE"
  9. 22  PRINT TAB(20);"OPEN MCNEMAR'S TEST (paired chi-square) OPEN"
  10. 25  PRINT TAB(20);"SCREEN";STRING$(36,205);"LOAD": PRINT
  11. 30  PRINT TAB(12);: INPUT "What is the name of the FACTOR to be tested?   ",F$
  12. 40  PRINT: PRINT "  ENTER the number of PAIRS in each category:"
  13. 50  PRINT: PRINT TAB(36);"CONTROLS"
  14. 60  PRINT TAB(28);"+ ";F$;TAB(42);"- ";F$
  15. 62  PRINT TAB(24);"VARPTRSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDBSAVESOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDCOLOR"
  16. 64  PRINT TAB(11);"+ ";F$;TAB(24);"CALL";TAB(39);"CALL";TAB(54);"CALL"
  17. 65  PRINT " CASES";TAB(24);"BLOADSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDBEEPSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUND<0xB4!>"
  18. 67  PRINT TAB(11);"- ";F$;TAB(24);"CALL";TAB(39);"CALL";TAB(54);"CALL"
  19. 70  PRINT TAB(24);"CLSSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDMOTORSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUNDSOUND'"
  20. 80  LOCATE 12,30: INPUT;"",AA: LOCATE 12,46: INPUT"",AB
  21. 90  LOCATE 14,30: INPUT;"",AC: LOCATE 14,46: INPUT"",AD
  22. 100  X=ABS(AB-AC)-1: X=X*X/(AB+AC): PRINT: PRINT
  23. 110  PRINT TAB(28);"CHI-SQUARE = ";X: PRINT
  24. 120  IF X>31 THEN P=0: GOTO 180
  25. 130  R=1.77245: S=1: I=1: K=SQR(X/2)*2/(EXP(X/2)*R): BD=3
  26. 150  I=I*X/BD: S=S+I: BD=BD+2: IF I>0 THEN 150
  27. 170  P=1-K*S
  28. 180  PLAY "O2 MB MS T120 L16D-P8L32GL16A-P8 L3 D-": PRINT TAB(19);: COLOR 0,7
  29. 190  PRINT TAB(32);"p = ";: IF P<9E-09 THEN PRINT "< 10 (-8)"; ELSE PRINT P;
  30. 200  PRINT TAB(60): COLOR 7,0: PRINT: PRINT: PRINT: PRINT TAB(10);
  31. 210  INPUT "Do you want to calculate another McNemar test? (Y or N)  ",A$
  32. 220  IF A$="y" OR A$="Y" THEN CLS: GOTO 20
  33. 230  END
  34.